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The generalized Helmholtz equations of relativistic multifluid plasmas can be integrated for 
axisymmetric equilibria in close analogy to the magnetic flux conservation law in ideal magnetohy- 
drodynamics (J. D. Bekenstein and E. Oron, Phys. Rev. D 18, 1809 (1978)). The results are, for 
each fluid component, two flux functions and a potential equation for the poloidal stream function. 
Ampere's equation for the four-potential A v is reduced to two coupled equations for the time-like 
and the toroidal component. So we have altogether four potential equations for a two-component 
plasma; they can be derived from a variational principle. 
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I. INTRODUCTION 



Observations of active galactic nuclei and other massive objects have increased the interest in equilibrium plasma 
models within a Schwarzschild or Kerr geometry. Even if the metric is given and all physical quantities are assumed to 
be independent of a toroidal angle ip and of time t, there is still some effort needed to reduce the whole set of equations 
for the electromagnetic field and the fluid quantities. This has been done within ideal magnetohydrodynamics (MHD) 
by several authors with various degrees of completeness and sophistication. A characteristic feature of the MHD 
description is the conservation of magnetic flux in a system co-moving with the center-of-mass four-velocity (of 
ions, essentially); this equation ("Ohm's law", Eq. (||) below) has been "integrated" for general stationary and 
axisymmetric systems by Bekenstein and Oron [Q who give also a basic discussion and some historical background of 
general-relativistic MHD. As a result, one can represent the electromagnetic field tensor completely by the particle 
flux nu 11 , the two Killing vectors associated with the translational symmetry in t and (p, and two "flux functions" 
which are constant along the poloidal stream lines (Eq. (|l^) below). Further reductions of the whole set of MHD 
equations and discussions of the astrophysical background have been given by Camenzind [g] , Mobarry and Lovelace 
, Nitta, Takahashi and Tomimatsu Q, and Beskin and Par'ev ||. Thus one arrives, as in the non-relativistic case, 
at a single potential equation for the magnetic flux function ^ (the covariant toroidal component A v of the vector 
potential A), together with some constraints. 

One problem with these MHD models is the large number of arbitrary flux functions - no dissipation mechanism 
has been included to reduce this number -, and a fluid picture may be questionable if the collision frequencies are too 
low. This, however, is not our point here, and we admit an ideal fluid picture on reasons of simplicity. The question 
is, however, whether the usual MHD theory is consistent for a really rotating black hole ||. In this case the metric is 
necessarily non-diagonal; the relevant element g tlp is the invariant scalar product of the two Killing vectors mentioned 
above 0, so it can by no means be transformed to zero. On the other hand, it couples the components At and A v in 
Ampere's law which will, in general, contradict the strong coupling of A t and A v in the MHD models (here they have 
to be functions of each other) . So we are led to the question of how to reduce the set of multifluid plasma equations 
and Maxwell's equations for stationary axisymmetric systems without the discrepancy mentioned above. 

In Section II we give a short account of the work of Bekenstein and Oron, leading to the strong coupling of A t and 
A v . The same solution, however, can be used for any generalized Helmholtz equation in ideal fluid descriptions. Two 
examples of these Helmholtz equations are derived in Sections III and IV, respectively: One refers to the ordinary non- 
relativistic MHD equation (momentum balance), the other to the relativistic multifluid plasma, assuming a constant 
temperature of all species. The latter example is used for the reduction problem in Section V. Here it is shown that 
the multifluid equations are reduced to one single potential equation for each species, and the poloidal components of 
Ampere's equation are integrated. A numerical solution of this full set of four potential equations would be facilitated 
by the existence of a variational principle; the solution could then be approximated with finite elements by minimizing 
the corresponding functional. So it may be interesting that such a functional exists, as is shown in Section VI. Finally 
we discuss the results in Section VII. 



II. FLUX CONSERVATION IN AXISYMMETRIC PLASMAS 



Let us assume a plasma configuration where all physical quantities, including the metric tensor components g^v, 
are independent of time t and of a toroidal angle (p. We choose a coordinate system (x^) with x° = ct (c = 1), 
x 1 = ip, and x 2 ,x 3 some poloidal coordinates. Assuming in addition that all physical quantities are invariant to the 
simultaneous inversion of t and ip - which is reasonable for any rotating equilibrium - the most general line element 
ds can be represented as follows || : 

(ds) 2 = g rs dx r dx s + g ab dx a dx b , (1) 

where the indices r, s run from to 1, and a, b from 2 to 3. We want to determine an electromagnetic field tensor 
Fpu which is consistent with Eq. (|l|), and which, in addition, obeys the condition of magnetic flux conservation in a 
medium moving with the Eulerian four- velocity u 11 . The latter condition is familiar from ideal magnetohydrodynamics 
(MHD) and means that a certain four- vector ("co-moving electric field") should vanish: 

= F MV u" = 0. (2) 

The field tensor F^, of course, should solve the homogeneous Maxwell equations; therefore it can be written as the 
curl of a four-potential A^ in the usual manner: 

Ffj,u = A y ^ A^ U: (3) 
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where is the partial derivative of the quantity in brackets with respect to x^. Finally we use the continuity 
equation 

(^Tjnun.^O, (4) 

with g the determinant of g^ v and n the particle number density in the local inertial rest frame. 

To solve Eq. (^) we remember that any skew-symmetric tensor can be represented by two four- vectors, 
and ( "co- moving magnetic field" ) [Q ; denoting all components in the local inertial rest frame by primes (with 
u ' = c = 1; u 1 ' = for i = 1, 2, 3) we define 

E' a = ; = 4 i= 1,2,3 

B 'o = ; B 'i = F jk' 

where in the last equation k) is a cyclic permutation of (1, 2, 3). In this coordinate system we use the Minkowski 
metric 

(s'nu) = diag(l, -1,-1,-1), 

and E %l — —E[. B 1 ' = —B^ are the local electric and magnetic fields, respectively. It is then easily seen that the 
covariant representation of E^ in the laboratory frame with any velocity field u^{x) is given by the left part of Eq. 
(0), while B^ is given by 



B^ = --£^u v F^, (5) 



where the totally antisymmetric Levi-Civita tensor e^uga includes a factor \J—g in order to be a tensor. Both vectors 
E^, B^ are obviously orthogonal to this corresponds to 2x3 independent components. Therefore we can represent 
the six independent elements of F M „ by these vector components, namely: 

F^v = E^u u - EyU^- e^ vea u B B a '. (6) 

The covariant equation (^) can easily be proved by writing it in the local inertial rest frame. Eq. (|^) is generally 
valid; in the case of flux conservation according to Eq. (^) it simplifies to 

F[iit e ftisga-u^ B . (7) 

Obviously the field tensor F M „ is then orthogonal not only to u M , as required by Eq. (||), but also to B^: 

B^F^ = 0. (8) 

The remaining task is now to construct B^ for axisymmetric equilibria. 

For this purpose we consider the two Killing vectors associated with these symmetries, namely: 

(#») ee (1,0,0,0), 
(f) = (m")= (0,1,0,0). 

In both cases they lead to a vanishing partial derivative of any physical quantity A in the direction of £^: 

= 0. (9) 

Let us assume that all components A v of the vector potential share this property (though a gauge transformation 
could destroy it). Then Eq. (||) leads to 

F rv ee ^F^ = -^A^ = -Ar lV , (10) 

where r = for £^ = fc M and r = 1 for £ M = m} 1 . The right-hand side of this expression is obviously non-zero only for 
v = a = 2 or 3. Differentiating Eq. (10) with respect to x b , b = 2 or 3, but b ^ a, we find 

Fra,b A ra fo F r b a . (H) 

To obtain the same expressions from Eq. ([?]) we write the factor yj—g of e^qo explicitely, with the remaining constant 
permutation symbol E^ vga : 
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Then our integrability condition from Eqs. (J?]) and ( |l~Lj ) reads as follows: 

= e rsab [(y/=g u s B») b - (V^9 u b B s ) b 

The right-hand side of this equation is understood with fixed and mutually different values for r, s, a and b. We can 
re-write it by restoring the summation convention with respect to the index a: 



(V-gu s B«) a -(V^u"B°) a = 0; 8 = 0,1. 
The general solution B 11 of Eq. ( [l2| ) with B^u^ = can be written as follows: 

B ll = au ll + b^ ■ a=-b^u,„ 



(12) 



where 6 M solves the same Eq. ([12]) as Since it is linear in we can solve it separately for the poloidal part b a 
and for 6 s . Ignoring now b a we can solve Eq. 
represent b s by a linear combination of the Kil 



for b s if the poloidal flow u a is not identically zero. It is useful to 
ing vectors, namely: 



(13) 
(14) 



where the coefficients C and (3 have to be determined suitably. The factor n has been included in order to take 
advantage of the mass conservation law, Eq. (^), where the replacement /i — > a is allowed due to the symmetries. 
Then we have a solution of Eq. (|l2| ) if C and /? are constant along the poloidal stream lines ( "flux functions" ) : 



u a C\ a = = u a f3, a , 

B" = -nC [K" - (K x u x )u"] , 



(15) 
(16) 
(17) 



The last equation has been obtained by inserting Eq. ( |l6| ) into Eq. ([?]). While C is an arbitrary flux function, j3 is 
fixed by the condition that F^ v is perpendicular to both u^,B^ or u^',K' 1 . Using Eqs. @, (§), ©, (||) and © 
we find 



= iW M „ = -(Aq^ + ZLV). 



(18) 



This equation can only be fulfilled if Aq and A\ are functions of each other and constant along the poloidal stream 
lines. This is indeed the case as can now easily be shown from Eq. (0) for /i = r and Eq. (fl0|): 



= F rv u" 



-A r _ a U . 



(19) 



This completes our particular solution for B M and F^ v if 6 a = 0. It coincides with the result of Ref. 01 (their A is our 
— (5). Why is a poloidal part of 6 M not possible? Our symmetry requires F rs = according to Eq. (|l0|); from Eq. (Q) 
we have 



T ^£T = -e rsab u a B b 



Here we see that for a non-zero poloidal flow it a the poloidal part of B^ must be proportional to u^, otherwise F rs 
would be non-zero. Therefore the solution as given in Eqs. (|l^) and (JT^) is unique, up to the specification of two flux 
functions, C and (3. 

Sometimes it is useful to re-write the result of this section in particular poloidal coordinates as defined by the 
poloidal stream lines, \& = const., and an angle-like coordinate 9 varying along the poloidal stream lines: 

x 2 = * ; x 3 = 6. 



From Eq. (19) we know that A r = Arj^Sf). so we may identify one of both components with itself, while the other 
component defines j3 according to Eq. (fill), e.g.: 
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Ai=\f ; /3 = -dA (y)/d*. (20) 
In this coordinate system we have from Eq. ( |l0| ) : 

-Pi 3 = ; F12 = — 1. 

Inserting here Eq. ( [L7| ) we find 

u 2 = ; C=f/(y^nu 3 ). (21) 

The remaining elements of F^ v in this coordinate system are then: 

F 0a = -0F la ■ F 23 ^((3u°~u 1 )/u 3 . 

Then our vector and field tensor elements F^ are, for a given geometry and flow field, completely determined by 
one single flux function (3. 

III. VORTICITIES IN MHD PLASMAS 

Flux conservation laws are generally expected if dissipative processes are negligible. Let us first discuss the 
non-relativistic case. The oldest example is the Kelvin/Hclmholtz theorem for a neutral fluid with pressure p = p(g), 
where g is the mass density. Euler's equation may then be written using the vorticity uj and the Bernoulli function 
U as follows: 

+LJXV + VU = 0, (22) 

and we obtain immediately Helmholtz's equation for the vorticity by taking the curl: 

-+Vx(uxd) = 0. (23) 

This equation is the prototype of any vector field u> — > f2 which is "frozen in", moving with the fluid velocity v, and 
which is the curl of another field, say V: 

dfl 

— + Vx Oxu = ; n = VxV. (24) 
ot 

To get an equation for V we integrate Eq. (]24|), introducing a scalar potential $: 

dV 

— + ftxv + V$ = 0. (25) 
ot 

For V — A and f2 = B we obtain 

E + vxB = 0, (26) 

where the electric field E is the usual expression in terms of the potentials A, $ (with c = I). Eq. ( |2^ ) is just 
the non-relativistic limit of Eq. (Q). For other ideal fluid models one may also find a Kelvin/Hclmholtz theorem, 
though the corresponding vectors V, l~i are more complicated. The momentum balance of ideal MHD theory refers to 
a "center-of-mass" fluid with total pressure p, total mass density g, center-of-mass velocity v, and the Lorentz force 
due to the electric current density j: 

dv „ . _ 
B ~dt + P = JX 

The picture of an idealized fluid requires not only zero resistivity according to Eq. (Pq) , but, more generally, zero 
entropy production. So we expect for ideal MHD theory a second Kelvin/Hclmholtz theorem and associated vectors 
V, fl; they seem to be unknown, but they can be constructed using appropriate Lagrangian and Lin variables p0[ . 
Here we use three Lin variables: the entropy per mass s and the two Euler potentials q\ (A = 1,2) of the magnetic 
field: 
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B = V<?i x Vq 2 - 



The constraints of entropy and magnetic flux conservation are then expressed by these three material invariants: 

ds dqi dq2 



0, 

at at at 



and the Ansatz for V reads as follows: 
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V-rVs-^vxVqx- (27) 



A=l 

The coefficients r, v\ are now determined in order to match Eq. (|25| ) with an arbitrary potential From Eq. (|2 
and its curl, the equation of O, we find after some rearrangements the following purely kinematical relation: 

dV dv dr ds 

— + ftxv = — +ivxv- — Vs + — Vr 

OT OT Cff Cff 

A A 

-vfg + x>t). <*> 

Here we insert our particular fluid model for dv/dt. The ideal MHD equation can be written as follows: 

dv f v 2 \ ?' 

— + u>x«=-V — + h}+ TVs + -xB, 
at \ 2 J g 

where h is the enthalpy per mass, T the temperature, and j the current density as determined from Ampere's law. 
The right-hand side of Eq. ( p8| ) is then just — V$, with 

i v 2 ds tt-^ dqx 

A 

provided the coefficients r, v\ obey the following equations of motion: 

dr „ dv^ 1 . dv2 1 . 

T =T ; -± = -j.Vq 2 ; -^ = ~J-V«i. 

dr dt g dt g 

So we find a second flux conservation law in ideal MHD; it refers to the following generalized vorticity: 

n = VxV = u - Vrx Vs - v "i x v< ?a- (29) 



The generalized Helmholtz equation ( |24| ) allows, of course, the "trivial" solution $7 = 0, corresponding to a potential 
flow for V , V = — VS 1 ; then Eq. (E7I) leads to the so-called Clebsch representation of v: 



-VS + rVs + Y v\Vq> 



A=l 

(For f2 we would need a further pair of Clebsch variables v\,q\ with dv\/dt — dq\/dt = 0). In contrast to the 
first Kelvin/Helmhotz theorem associated with Eq. (|2^) we have no simple advantage from this second vorticity law; 
in particular, the equations of motion for r and v\ have to be solved, in addition to the conservation laws for s and 
q\ and (not shown here) the Bernoulli equation for S. The situation, however, becomes more transparent if we leave 
the MHD description and treat the plasma as a fluid of different species j (j — e : electrons, j = i : ions of any kind) 
interacting only via the electromagnetic field, with the following momentum balance for each species: 

Qi (§1 + v i ■ v ) v i + ^Pj = e i n 3 ( E + v i x B ) 



G 



There is an obvious formal bridge to the MHD description for a two-component plasma: Neglecting the electron mass 
(m e ~ g e — > 0) we obtain Vi as the center-of-mass velocity v; assuming, in addition, quasi-neutrality with simply 
charged ions (n e = rii = n) we can replace v e as follows: 

3 

v e = v- —. 

\e\n 

Adding both momentum equations gives then the correct ideal MHD balance with p = p e + Pi, while the momentum 
equation of electrons leads to a generalized Ohm's law replacing Eq. (p6|): 

E + vxB= —(jxB-Vp e ). 

\e\n 

This "modified" MHD theory is the bridge mentioned above, but a simpler form of the two Kelvin/Hclmholtz theorems 
is obtained if we come back to the multfluid description, now within general relativity. 



IV. VORTICITIES IN ISOTHERMAL PLASMAS 



The flux conservation theorems associated with perfect fluid models depend crucially on the equation of state 
for the pressure, e.g., p = p(g,s). The non-relativistic case in standard textbooks on neutral fluids usually assumes 
constant entropy s throughout the fluid volume, leading to the usual Hclmholtz equation, Eq. (^3|). The corresponding 
equation in general relativity has been derived by Taub |l3| . If the entropy of a non-relativistic neutral fluid varies 
in space, the situation is more subtle; it has been discussed already by Eckart |L4j and later in Ref. |J. The result is 
that the vector V whose vorticity flux is conserved differs from v by (— rVs), the second term on the right-hand side 



of Eq. (27). For a non-relativistic multifluid plasma Eq. (|27j) is replaced by 



leading to a "three-circulation theorem" corresponding to the three constituents of V (A is the vector potential). In 
general relativity we have, instead of Eq. ( p5| ) (from which the general Helmholtz equation (|24|) follows), the equation 
corresponding to Eqs. (j|) and (||) for each species (suppressing now the species index j): 

uTl^ = 0, 

The vector V v for a multifluid plasma with s — const, has been given by Lichnerowicz (jl5| and Carter |L6| ("single 
constituent perfect fluid", his example (c) in §4), and for varying s by Ref. |Tl]] , namely: 

Vu = —Ua + rs t u H A^, 

g m 

where a is the relativistic enthalpy per volume, and the covariant Eulerian four- velocity (with c = 1). One of 
Carter's results refers also to a neutral fluid with two constituents and varying entropy, but the resulting canonical 
momentum per volume (= gV^) reads, in our notation, as follows (from Eqs. (4.25), (4.26), (4.32) and (4.33) of Ref. 
@): 

gV^ = 71717, = f7U M' 

so the term rs^ is missing. 

Here we consider a different physical situation which may be of astrophysical interest: We assume that a radiation 
field acts like a heat reservoir for electrons and protons, maintaining a constant temperature of them. Then the term 
rs iAI is again absent, and we find the relativistic Helmholtz equation for the ordinary canonical vorticity of both 
species. In this case, a will turn out to be the free relativistic enthalpy per volume. We start with the material 
energy-momentum-stress tensor T M " of an ideal electron or ion fluid without denoting the species index explicitely. 
Since both fluids are only coupled via the electromagnetic field by the Coulomb/Lorentz force, we can write the 
energy-momentum balance for both species as follows: 

T„ = en F^u". (30) 
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The semi-colon indicates the covariant derivative. The coupling of both species by the right-hand side of Eq. ( J30| ) 
implies now that the magnetic flux is not conserved, neither for the electron nor for the ion fluid. The tensor " for 
an ideal fluid is well-known (see, e.g., Ref. |l2[): 

The last term in this equation is the scalar pressure p in the local inertial rest frame times the Kronecker symbol, and 
the scalar a depends on the equation of state. Using also the mass conservation, Eq. ([!]), we obtain 

T M u - u = Qu v (jj U i^J ~V,n- 

This four- vector must be orthogonal to as is also the right-hand side of Eq. (|30|). With the normalization of 
u M {u^u^ — 1) we find then 

= v , v = gu v Q - u"p, v , 

d(-) = -dp, 

where the differential d means variation along the path of a fluid element. In the co-moving inertial rest frame we use 
the Gibbs-Duhcm relation for the free enthalpy p, per mass: 

dfx = — s dT H — dp. (31) 
Q 

Ignoring temperature variations, and including the relativistic rest energy gc 2 (with c ^ 1 for the moment) the 
resulting expression for a is then as follows: 

* = Q (l + ^). (32) 

To derive a flux conservation law we insert the vector potential for F^ v in Eq. ( |30| ) and put all terms to the left-hand 
side: 

a \ e e 1 

g J m m g 

To eliminate the pressure term we calculate the partial derivative from Eqs. ([}]]) and J32] ) with dT — (and this time 
c=l): 

where in the last step we used again the normalization of u v . This is now just the term leading to flux conservation 
for the canonical vorticity of each species; we define 

V M = -Un + —An, (33) 
q m 

% v = V Vili - V KV = V v ,ix - V^, (34) 

and we find 

% v u" = 0. (35) 

The solution of Eq. ( |3^ ) for axisymmetric equilibria is now simply obtained from the previous section: We replace 
there F^ v by and by V^. The mean velocities of electrons and ions, however, are usually different; the fluxes 
of their general vorticities £1^ are therefore conserved in different frames. It is interesting to rewrite Eq. (|35|) for the 
spatial components of the canonical velocity, V 1 , in the special-relativistic case (no gravity), namely: 

dV 

— + tlxv + VU = 0. 

at 

where U{= cVq) is the relativistic Bernoulli function, and (= — £ljk) are the spatial vector components associated 
with O^y. This equation is now again of the same form as Eq. (p5|), and we recover Eq. (^4|) as the prototype of any 
special-relativistic generalized Helmholtz equation. 
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V. MULTIFLUID PLASMA EQUATIONS FOR AXISYMMETRIC EQUILIBRIA 

Let us use a general poloidal coordinate system to solve the mass conservation law, Eq. (0), for each species 
separately by introducing an appropriate stream function Xi namely (e lab is again the permutation symbol, here for 
spatial indices): 

" - ' e lab X , b . (36) 



y/— g n 

The poloidal stream lines are then given by x — const., and Helmholtz' equation ( |35| ) in the symmetry plane, 
(X = r = 0, 1, reads as follows: 

Q = -Vr,„u v = -V r , a u a . 

Similarly as in Eq. ( |l9| ) we conclude that the components V r can be any flux functions with respect to x : 

Vr = V r ( X ). 

These two flux functions and the corresponding components of A v fix two components of u u for each species up to a 
factor g/a, namely: 

5 f Vr(x) ~ -A r ) . (37) 



a \ m 

Assuming that x an( i A v are given elsewhere we may read the normalization condition for u v as an equation for n: 

1 = g rs u r u s + g ab u a u b . (38) 
So all fluid quantities besides x are determined by Eqs. ([36]) - fl38|), and we find all elements of f2 M „ except fl a b- 

Q rs = ■ fl ra = -V r , a = -V^a, (39) 

where the prime means differentiation with respect to x- The general solution for f^j,, however, can be obtained from 
Eqs. ( |T7j ) and (18) with appropriate changes of notation, namely: 

= nC{ X )e^ ea u e K°, (40) 
= -V6(x)/Vl(x). (41) 



Comparing this with the results above (Eqs. (g6|) and (|39| ) ) we find the flux function C(x): 

C(X) = -Vl(x). (42) 

The equation for the stream function x is then obtained from Eq. ( ^p[ ) with /.i = 2 and v = 3, where the left-hand 
side is determined according to the definitions (HJ) and (^3|); the result is then the following: 



"23 = V=g nV;( X )u r , (43) 

n 23 =(-u 3 ) - ( -u 2 ) +-F 23 . (44) 
\Q J , 2 \Q J ,3 m 

The complete set of fluid equations (|3^) - (^8|) and (^3|) - ( f4~I| ) for the unknown variables u r , n and x is written for 
any poloidal coordinate system, thus allowing any number of particle species. To solve finally Ampere's equation in 
the poloidal plane we are then free to use particular coordinates and a particular gauge of A v . To be consistent with 
the usual notation we denote the projections of the -lines onto the poloidal plane the lines ^ = const., where the 
flux function ^> may be identified with a "radial" coordinate x 2 as in Sec. II, and the stream function ~ / of j a is 
denoted as a flux function, / = I^). The continuity equation for in the poloidal plane is then solved as follows: 

f=0 ; 4ttV^ f =/'(*), (45) 

where the prime of I means differentiation with respect to 'J. The flux function I(JSf) is, of course, not independent 
from the stream functions x of the different species. From 
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E 



enu' 



and Eq. (45) we find the following relation: 



E 



ex 



47T 



const. 



(46) 



where the sum over j is the sum with respect to the different particle species. A convenient gauge of A v is, as in the 
non-relativistic case, the condition that A is tangential to the surfaces ^^const. of the current lines: 

A W = 0. 

This equation can usually be fulfilled by an appropriate gauge function since \& is a radial coordinate; for \& — x 2 it 
reads 



A 2 = or F 23 = A 3>2 . 
Let us now start with Ampere's equation for A v : 

y—g g^g v ° (A a , e - A e . a )] v = -An^d f 



(47) 



(48) 



In the symmetry plane, fi = r = 0, 1, this equation is decoupled from the poloidal components of A v \ inserting the 
fluid quantities for j r we have then the following set of two equations for the components A s , s = 0, 1: 



•//"•''.!-,; „ = ^v~g g rs E e n - ( v > M - - 

' * — ' a \ m 



(49) 



In the poloidal plane, Eq. ( f48| ) can be re-written with indices a, 6, c, d which run from 2 to 3 only: 

[^dg ab g cd F bd } c = -An^f 

The left-hand side is easily evaluated due to the antisymmetry of F^', we introduce the determinants of g a b and g r 
explicitely: 

g P oi = det(g ab ) ; g sym = det(g rs ), 
then we obtain the following form of Ampere's equation in the poloidal plane for a = 2: 



,9.s 



9poi 



' J? 



2.'! 



and for a = 3: 



gsyr, 
9pol 



■F. 



2:i 



-£-K^f=gj< 



-Ait^f^g j . 



Inserting here Eq. (p5[) we identify the expression in brackets as a flux function, namely: 



9syr< 
9pol 



■F 23 = /(*). 



(50) 



(51) 



(52) 



In the non-relativistic case this is simply the covariant toroidal component of the magnetic field. Finally, the component 
A 3 is determined from Eqs. (p2) and ([47|); it is not a flux function, and it is not needed in the remaining set of equations. 
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VI. A VARIATIONAL PRINCIPLE 



The numerical solution of our potential equations for x, Eqs. ( B ) - (f44|), and A ri Eq. (49), is simplified by the 
existence of a functional of x and A r which is stationary in equilibrium. We need, however, still another quantity 
whose variation leads to the normalization condition, Eq. ([58]). So we look for a functional W of three quantities, 
W — W{a, Xi Ar) sa y, whose independent variations lead to the equilibrium conditions. The particle density n (which 
is not varied) , and a of each species are then obtained afterwards from a combination of the normalization condition 
an an equation of state according to Eq. J3^ ) . 

Let us start with the normalization condition, where u a and u r are given by Eqs. (|3^) and (j37[), respectively. 
The latter equation gives 



1 Q CTl 

5u r — u r 5a H — Vr(x)$X 5A ri 



(53) 



while 5u a depends on Sx only: 



Su a = ^^e lab (S X ).b. 
■J—g n 



It is then easily realized that W could be of the following form: 



(54) 



s u r u s — 1) + 



where the terms indicated by dots are independent of a, and the integration is done in a fixed region of the poloidal 
plane (x 2 , x 3 ). It is a remarkable effect that in the expression above the usual invariant u v u v is replaced by u a u a — u r u r 
which is only invariant under transformations in the poloidal plane. A similar replacement will be needed in the 
invariant which produces Maxwell's equations in vacuum. Using our solution for F23, Eq. (p2[), and F rs = 0,F ar — 
A r a , we find 



1 

167T 



1 7 2 (*) 1 

V = —zz , r - ~ ^zg rs g ah A r , a A s , b . 



8vr (-g sym ) 87T' 

Changing now the sign of the last expression above, we are led to the following functional: 



W(a, X ,A r ) 



d X\J — g 



X] -j^9abU a u h - g rs u r u s - 1) 
3 

— 7 — — H g g A r a A s b 

871" (~gsym) 87T 



(55) 



Variations with respect to x an d A r are now done by eliminating derivatives of 8x and SA r by partial integrations, 
assuming that Sx and SA r vanish at the boundary. Furthermore, we have to vary /('J) according to Eq. (E6h , namely: 



= (-4we)6x- 

The total variation of W(a,x, A r ) is then obtained with the following result: 



(56) 



SW 



J d 2 xV^ff|y^ ^(g a bU a u b + g rs u r u s - l)5a 



(_ 1 n fee + e/(*) + g sym gu r v;(x) ) S X 



(V=99 rs 9 ab A Sib ) ta 



SA r 



(57) 
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Ax = /--^ y2(g ab x,b - gbbX,a),a, (58) 
a 

g a b = , — 2 ffafe' (59) 

The vanishing factor of 5A r is easily identified with Ampere's equation in the symmetry plane, Eq. ([l9|). To identify 
the vanishing factor of 5\ we note that from Eqs. m4), (p3), and (pfl) we have 



mn 23 = ^2(gabX,b ~ gbbX,a),a + ej — 



0,6 v 5sym 

Inserting this in Eq. ( p3| ) with a slight rearrangement, we find 

A X + e/(*) + g sym gu r v;(x) = 0. (60) 

This is just the condition that W is stationary with respect to variation of x, assuming that g sym is finite at this 
point. 



VII. DISCUSSION 



A plasma equilibrium near a rotating black hole has been considered in the ideal fluid picture. The usual MHD 
equations imply two flux conservation laws. One is the well-known conservation law of magnetic flux (Ohm's law with 
vanishing resistivity); it leads to two constants of motion for stationary axisymmetric systems: The covariant time-like 
and toroidal components, A t and A v , of the vector potential are constant on the poloidal stream lines of the plasma 
bulk velocity. This well-known fact is in contradiction with the coupling of A t and A v for a metric with g tlfi ^ 0, as is 
appropriate for a rotating black hole. Simple Grad-Shafranov type MHD equilibria (see, e. g., Ref. jl7| are then ruled 
out in this case. The second flux conservation law of MHD has been derived in Section III for the non-relativistic case 
(Eqs. ( p^ ) and ( p7| ) - (|2^)); it is, however, not as simple as the magnetic flux conservation law. A more reasonable 
description for g tip ^ is given by the multifluid equations: They can be cast into the form of Hclmholtz equations 
for each fluid component (jll); they are particularly simple for isothermal plasmas, as is shown here (Eqs. ([53|) - (|35|)). 
Since the form of these equations is exactly the same as the magnetic flux conservation law, we obtain two constants 
of motion for each species of an axisymmetric equilibrium, the time-like and toroidal components of the canonical 
velocity V v . The component V t is the relativistic Bernoulli function, and V v is the canonical angular momentum 
per mass of a fluid particle of a particular species; this has to be expected for axisymmetric equilibria if the fluid 
components interact only through the electromagnetic field. The poloidal components of Ampere's equation can be 
integrated, too, by adjusting the coordinates to the lines of the electric current, \& = const.; the relevant flux function 
1(4"), Eq. (|52|), corresponding to the toroidal magnetic field, is now simply related to the stream functions x °f the 
different plasma components according to Eq. (|46|). The whole set of potential equations can now be summarized 
in a functional W(a, x, A r ), where a is the free enthalpy density of a particular species, x the stream function of its 
poloidal velocity, and A r stands for A t and A v . The total variation of W produces then the nomalization condition 
for the Eulerian four- velocity of each species, the potential equation for x (Eq. (|6(i|)), and Ampere's equations for A t 
and A v . 

It is interesting to consider possible solutions of these equations in a given geometry with g tip ^ like Kerr's 
metric. Equilibria with poloidal velocity fields arc of interest from an observational point of view, because they are 
able to exchange mass, angular momentum etc. between inner and outer parts. Plasma models with pure poloidal 
velocity fields, however, are not possible; the reason is that the toroidal velocity u v is proportional to (V v — (e/m)A v> ) 
(Eq. (|37|)), where V v is constant on the stream lines of the particular species which is considered, but A v generally 
not. Equilibria with pure toroidal velocities are possible for an electron-positron plasma. In this case Eq. j60[ ) is 
solved trivially with x = 0, V r = const., = 0, and the velocity components u t , u v become equal but opposite 
in sign, up to a constant V r ; they are determined from Ampere's equation (^]) which is highly nonlinear due to 
the normalization condition. This solution, however, seems to be artificial, and an acceptable solution will exhibit 
inevitably also poloidal velocity components. 
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